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Abstract. An approximate approach for the description of the two-photon 
bremsstrahlung emitted by a relativistic projectile scattered in a spherically-symmetric 
field is developed. It based on the accurate treatment of the analytical structure of 
the singularities in relativistic one-photon free-free matrix elements. For the case 
of a Coulomb field the analytical expressions for the amplitude and cross section 
of process are presented. Numerical results obtained within the framework of the 
proposed approach are compared with the available experimental data and with the 
results of simpler theories. 



PACS numbers: 31.15Md, 32.80Wr 
1. Introduction 

In this paper new results from the relativistic theory of two-photon bremsstrahlung 
(2BrS) of a projectile electron scattered in a spherically-symmetric static field are 
reported. We formulate the approximation suitable for effective analytical and numerical 
treatment of a rather complex relativistic free- free two-photon matrix element. The 
developed formalism is applied to the 2BrS process in a point-Coulomb field. Results 
of numerical calculation of the 2BrS cross section for several geometries of the emission, 
incident electron energies and photon frequencies are presented. 

The approach, which is described in this paper, based on the use of a so-called 
'delta'-approximation. This method, initially formulated for a non-relativistic 2BrS [1]- 
[4] and extended later to the case of relativistic collisions [5], is based on the assumption 
that to a great extent the behaviour of the free-free two-photon amplitude is defined 
by the contributions of the delta-singular parts of the one-photon matrix elements from 
which the 2BrS amplitude is constructed. 

Transitions of such type are encountered not only in the 2BrS problem but 
in a number of other physical processes in which the collision of a projectile with 
a target is accompanied by emission/absorption of photons (see, e.g., reviews [6- 
8]). The list of such processes includes, in particular, (a) the bremsstrahlung- type 
phenomena in an external field (many-photon spontaneous bremsstrahlung, laser- 
induced bremsstrahlung), (b) various inelastic processes, when the emission/absorption 
of the photon is accompanied by simultaneous excitation or ionization of the target, 
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(c) Compton scattering from many-electron atoms, (d) many-photon ionization of atoms 
and ions. 

In the case when perturbation theory in photon-projectile (or photon-atom) 
interaction is used to analyze the above-mentioned processes, the corresponding 
amplitude M. can be represented in terms of the compound matrix element which 
contains the radiative free-free matrix element between the intermediate virtual states 
and in which the integration over the intermediate momentum is carried out. The 
general form of such matrix element is 

• M = E/ dp b ^ b ( Pb )^(k,e). (1) 

Mb 

Here D{^^ (k, e) stands for the matrix element of the one-photon free-free transition 
between two relativistic states of a continuous spectrum with momenta p a and Pb and 
polarizations /i a and /ly,: 

D^(k,e) = J dr*£>t(r)(ea) exp (-ikr) tf£>(r) . ( 2 ) 

Here ty„ (r) are the bispinor wavefunctions corresponding to the out- (the upper index 
'+') and to the in (' — ') scattering states, the symbol f denotes the hermitian conjugation, 
ex. = 7°7 where 7 , 7 are the Dirac matrices. The vectors hk and e denote the photon 
momentum and polarization. We choose the gauge in which ek = 0. The CGS system 
is used throughout the paper. 

The factor ^^(pb) represents the contribution of all the remaining processes 
related to the collision. In particular, in the case of 2BrS process, ^^(pb) contains 
another free-free matrix element (see general expression (jUJ) for the 2BrS amplitude). 

Relation (J2J) defines four matrix elements the type of which depends on the 
asymptotic behaviour of the wavefunctions of the initial and final states. In the vicinity 
of the points defined by the relations 

Pa-Pb-ftk = 0, p & - |p b + hk\ = , I p a - hk\ - p h = , (3) 

the matrix elements D^t^ (k, e) contain the singular terms of the two essentially different 
types [5]. Firstly, there are well known pole-like singularities which are responsible, in 
particular, for the infrared divergency of the one-photon bremsstrahlung amplitude (see, 
e.g., [9-11]. In this case both momenta, p a and pb, lie on the mass surface, so that the 
exact equalities (JHJ) are inconsistent with the energy conservation law, which can be 
written as e a = + hu (here are the energies of the projectile in the states 'a' and 
'b', huj is the photon energy). Hence, in this case the equality signs must be understood 
as the limits — > 0. 

The singularities of the second type (we call them 'delta' singularities) contain terms 
proportional to the delta-functions 5(p a — pb — ftk), 5{jp & — |Pb + ^k|) and <5(|p a — hk\ —ph) 
(note, that the first of these is the three-dimensional 5-function, whereas the two others 
are one-dimensional). In the compound matrix elements (0) one of the states 'a' or 'b' 
(or both, if the process includes the emission of two and more photons, or if the photon 
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emission occurs between two virtual states) describes the off-mass surface particle. 
Therefore, the 'delta' terms add a non-zero contribution to the amplitude (JIJ. 

The approach for the approximate treatment (the 'delta'-approximation) of the 
compound matrix elements of the type (JTJ) is based on the assumption that to a great 
extent the behaviour of M. is defined by the contributions of the delta-singular parts of 
the one-photon free-free matrix elements D^t^ (k, e). Initially [1-3] it was formulated 
in connection with the problem of two-photon dipole free-free transitions of a non- 
relativistic projectile moving in an external field of atomic target. In the cited papers 
as well as in [4, 12] it was demonstrated that over wide regions of the incident electron 
energies and the photon frequencies the approximate approach produces nearly the same 
results as the rigorous treatment for both dipole [13-16] and non-dipole [4,17] photons. 

Within the framework of the non-relativistic dipole-photon approximation it was 
demonstrated that for the correct evaluation of it is necessary to take into account all 
singular terms of the free-free matrix elements, and in many cases the contribution of the 
delta-terms to the integral on the right-hand side is highly noticeable, if not dominant. 
Such an analysis has been carried out in connection with many-photon ionization 
(detachment) [18-21], the single-photon bremsstrahlung process with simultaneous 
ionization of the target [22], the process of two-photon bremsstrahlung [23,24], and 
the process of spontaneous bremsstrahlung in presence of a laser field [25]. In non- 
relativistic non-dipole approximation the analogous treatment has been applied so far 
to the free-free compound matrix element which characterizes the amplitude of two- 
photon bremsstrahlung [3,4]. 

A detailed review of the achievements of the non-relativistic theory of 2BrS can be 
found in [4,5]. Here we just mention the main results obtained in this field. 

Theoretical activity was stimulated by a series of experiments [26-33] in which the 
data on the 2BrS cross section and angular distribution were obtained for relatively 
high energies of the incident electron and emitted photon (up to 70 keV). Apart from 
the early work of Smirnov [34], where the 2BrS process was treated in terms of the 
relativistic first Born approximation (see also [35]), theoretical investigations were based 
on second order non-relativistic perturbation theory. Within the framework of this 
approximation the exact analytical formulae were obtained for the compound free-free 
matrix element in the Coulomb field. It was done in the dipole-photon case [13-15,23,36] 
and with the retardation effects included [4, 17] as well. For neutral atomic targets the 
calculations of the 2BrS spectrum were performed within the frame of the potential 
('ordinary') bremsstrahlung model [1,12,37] as well as with the 2BrS emission due to 
the 'polarizational' bremsstrahlung mechanism taken into account [38-40]. 

In connection with the experiments [29-32], where the 2BrS process was studied 
for 70 keV electrons scattered by various many-electron atoms, it was noted in several 
publications [1,3,4,16,17,41] that the role of the retardation effects in the 2BrS process is 
much higher than in the conventional one-photon bremsstrahlung. In these papers it was 
also demonstrated that that for the energies of the electron and the photons within the 
tens of keV range (i.e. as in the cited experiments) it is the retardation and relativistic 
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effects which are mostly important for the formation of the 2BrS spectrum rather than 
the effect of screening due to atomic electrons. The inclusion of the latter effect in 
the calculation scheme mofifies the result on the level of several per cent (see, e.g., [1]) 
whereas the account for the radiation retardation increases the 2BrS cross section by 
the order of magnitude [4, 17]. Therefore, to obtain a reliable theoretical result one has 
not only go beyond the frame of the non-relativistic dipole-photon approximation but 
to consider the fully relativistic treatment of the problem. 

The relativistic treatment of the compound many-photon transitions is much more 
complicated than its non-relativistic analogue from both analytical and computational 
viewpoints [42-53] . In the cited papers the theory and the numerical results for various 
processes described by the two-photon bound-bound and bound-free transitions are 
presented. In contrast, the exact relativistic treatment of the free-free two-photon 
transitions has not been suggested so far. Till recently all theoretical considerations 
did not go beyond the framework of the plane wave first Born approximation [34, 35] 
and the soft-photon limit [54]. In the recent paper [5] the relativistic formalism of 
'delta '-approximation for treatment of many-photon free-free transitions was developed. 
Its application was illustrated by the evaluation of the delta-amplitude of 2BrS 
process. Excluding the papers [34,35,54] there were no numerical results of relativistic 
calculations of 2BrS cross sections. 

In order to fill this gap we extend the approximate treatment of the relativistic 
2BrS developed in [5] to the case of the point Coulomb field, —Ze/r, and perform a 
numerical analysis of the role of retardation and relativistic effects in the formation of 
spectral and angular distributions of 2BrS for the conditions of the experiments [29-32]. 
This is done within the framework of the 'delta'-approximation, described in more detail 
in [5] and in section |21 below. Because of the absence of the closed analytic expression 
for the relativistic scattering states wavefunctions we used the Furry-Sommerfeld-Maue 
(FSM) wavefunctions (e.g. [9]). The applicability of FSM approximation in view of 
available experimental data is discussed in section l2~3l 

In spite of the fact that the validity of the 'delta'-approximation has not been 
formally proved there is a number of indirect confirmations of the applicability of 
this approach. In particular it was demonstrated (both in the non-relativistic [3] and 
relativistic [5] cases) that the general expression for the 'delta'-amplitude of the 2BrS 
process correctly reproduces the important limiting cases: the plane-wave first Born 
approximation and the soft-photon approximation. Also the 'delta'-approximation, 
being applied for the description of a number of mentioned above processes, yields 
reasonable numerical results against the results of the exact calculations. In the dipole- 
photon regime it was done for the two-photon bremsstrahlung [1,2,12,24], two- and 
three-photon detachment of electrons from negative ions [55]. Beyond the dipole-photon 
approximation the comparison was carried out for the 2BrS process of a non-relativistic 
electron [3,4]. 

This approach, although being approximate, allows us to evaluate effectively the 
principal parts of the free-free two-photon matrix element in the relativistic domain 
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Figure 1. Diagrammatical representation of the 2BrS process. The solid lines 
correspond to the relativistic projectile moving in an external field of the target (the 
Furry picture, e.g. [9]). The transition from the initial state '1' into the final state '2' 
via the intermediate (virtual) state 'n' is accompanied by the emission of two photons 
7i,2 (the dashed lines). 



with much less analytical and computational efforts. The method can be generalized to 
embrace the n-photon (n > 2) radiative free-free transitions. 

The formalism is described in section El In section |3] we present the results of 
numerical calculations and carry out the comparison with available experimental data 
and the results of simpler theories. 



2. Formalism 



2.1. 2BrS amplitude and cross section: general formulae 

The 2BrS process during the potential scattering of an electron (the charge is labeled as 
e and the mass as m p ) is a transition of the projectile from the initial state V\ = (p 1; fii) 
to the final state v 2 = (p 2 , fi 2 ) accompanied by the emission of two photons. Each of 
the photons (j = 1,2) is characterized by the energy, hujj, momentum ftkj and the 
polarizational vector e^. The energies of all particles satisfy the energy conservation law 

si = e 2 + hujx + kuj 2 • (4) 



where £\p = ypl 2 c2 + ?^pC 4 are the particle energies in the initial and the final states. 

The four-fold cross section differential with respect to u>i i2 and to the solid angles 
of the photons emission f2 12 is given by 

dV _ a 2 eie 2 p 2 uiu 2 \- \- f . |2 / \ 

do^d^dfla "2 Pl (2vr)6 2-, J^™ 1 ** ^ ' [b) 

Ai,A 2 A»i,M2 



Here a is the fine structure constant, the integration is performed over the solid angle 
of the scattered particle, the summations are carried out over the polarizations of the 
photons (Ai j2 ) and of the projectile (^1,2) • 

The total amplitude of the process, A4, is described by two Feynman diagrams as 
presented in figure ^ Each diagram corresponds to the compound matrix element of 
the two-photon transition from the initial state '1' via the intermediate (virtual) state 
£ n' to the final state '2': 

ka \ a 1 1 KA sr^ (p 2 Ai 2 |e 2 7e" ik2r |n)(n|ei7e~ ikir |pi / ui) 
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Here (p2A^| = ^LmI ( r ) 7° an d |piA*i) = ^piMi( r ) are the projectile's initial- and final- 
state wavefunctions. The sum is carried out over the complete set of states \n) of the 
Hamiltonian H = cap + 7°mc 2 + V(r), and includes the contributions of the positive- 
energy (e n > 0) and the negative-energy (e n < 0) states, the index n includes all 
quantum numbers which characterize the intermediate state. The term .Mpi] can be 
obtained from M.[\2) (the ordering of the subscripts inside the square brackets indicates 
which of the photons was emitted first) by exchanging (u>i, ki, e x ) «-> (dj 2 , k 2 , 62). 

The bispinor wavefunctions \?p^(r), which describe the scattering states with the 
asymptotic momenta p, can be expanded in the partial-wave series over the functions 
$? e ji m (r) characterized by energy, e, total angular momentum, j, orbital momentum, I, 
and projection of the total momentum, m [10]: 

= ^ E (4>p) ^( n p)) e± ^ i(£) *^ r ) • ( 7 ) 

jlm 

Here the notation n a stands for the unit vector along the direction a, Qji m (n) denotes a 
spherical spinor defined as in [56], is a unit two-component spinor corresponding to 
the spin projection \i. The quantities 6ji(e) are the relativistic scattering phaseshifts. 

The non-relativistic analogue of Mu 2 ] was investigated analytically (in the case of 
a point Coulomb field) [4, 13-15, 17, 23, 24, 36] and numerically (the cited papers and 
also [1,12,37,57]). 

Contrary to the non-relativistic case, the evaluation of the right-hand side of A^[i2] 
from (jnj) can be envisaged only by means of numerical calculations. Even in the case of 
a point Coulomb field one can hardly anticipate any real progress in applying analytical 
methods. This is due to both the complexity of the analytic structure of the relativistic 
Coulomb Green's function (e.g., see the discussion in a recent review [6]) and to the 
absence of the closed analytic expressions for the scattering states wavefunctions (J7J). 

The direct numerical computation is also a challenging problem. So far real progress 
has been achieved only in the fully relativistic treatment of the bound-bound and bound- 
free compound matrix elements (see, e.g., [49-53] and references therein). Even though 
it is possible, using (JJJ) and the corresponding expansions for the Green's function and for 
the operators ei !2 7e _lkl ' 2r , to write down the partial wave expansion for the amplitude 
A4, it is still a formidable task of computing the obtained multi-fold series. 

Alternatively, to compute the exact relativistic 2BrS amplitude one can apply the 
method due to Sternheimer [58] and Dalgarno and Lewis [59]. In this case, instead 
of the direct evaluation of the sum over the intermediate state n, one can evaluate 
.M[i2] as a matrix element of the operator e 2 7e _lk2r between the final state {P2M2I and 
some auxiliary function which is the solution of properly defined inhomogeneous Dirac 
equation. In the relativistic domain this method was successfully applied to the bound- 
free two-photon transitions [49], and to the bound-bound ones [52,53]. On the basis of 
preliminary analytical and numerical work which we have done (but do not report on it 
in this paper) it seems feasible to compute the 2BrS amplitude by this method. 
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However, the aim of this paper is to apply another approach, which, although being 
approximate, allows one to carry out the analysis of the characteristics of the relativistic 
2BrS with much less analytical and computational efforts. For the sake of convenience, 
in the next section we briefly outline the idea and the basic formulae of the approximate 
method [5] which is used further in the paper for the calculation of the 2BrS cross section 
of a relativistic electron in a point Coulomb field. 



2.2. 2BrS amplitude in the 'delta' -approximation 

To introduce the 'delta'-approximation let us start with separating the contributions of 
the positive- and negative-energy parts of the electronic propagator to the total matrix 
element M.\\2] from (JOJ): 

A* [12] = [M [12] ] + + [M [12] ]_ (8) 

To avoid the unnecessary complications we assume that the positive-energy 
spectrum of the electron in the external field V(r) does not contain the bound states, 
so that all intermediate states with e n > belong to the continuous spectrum. Hence, 
the sum on the right-hand side of (jBJ) can be understood as the sum over the bispinor 
polarizations fi and the integral over the momenta p. Hence, the amplitude [.M[i2]] + 
can be written as follows: 

L*W + -]LJ (2^3 e „ + /i Wl - ei -io - ^ [12] + 1/[12] • (9) 

The term J[i2] (which appears if one uses the standard rule to extract the imaginary 
part of the integrand's denominator) describes a two-step emission process in which the 
energy of the intermediate state satisfies the energy conservation law: e v = e\ — hw\. 
This quantity is given by the following integral 



J [12] 



pe 



8n 2 h 3 c 2 



E / dfi P ^ ±) (k 2 ,e 2 )^ 1 + )(k 1 ,e i ; 



• (10) 

=£1 — ?KJ1 

We note that because of the angular integral over Q p the result of integration in (JTUJ), 
as well as in Q, does not depend on the type of basic set (i.e. '+' or ' — ') chosen to 
describe the intermediate continuum states ^(r). 

The approximation, briefly described below (for the details see [3,5]), concerns the 
evaluation of the first term, R\i2], from Q which reads 

where v. p. indicates the principal value integration with respect to the pole coming 
from the integrand's denominator. The free-free matrix elements Di^vifai&z) an d 
^ (ki,ei) correspond to the radiative virtual transition. The momentum p of the 
intermediate state v is not fixed by any conservation law. The delta-singular terms 
in Di^^\k. 27 e 2 ) appear (cf. eq. Q) if p satisfies either | p 2 + ^k 2 1 — p = or 
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|p — Kk 2 \ — P2 = 0, which define two spheres in the momentum space. In the vicinity 
of these spheres the matrix element has pole-like singularity which must be treated in a 
principal-value sense when carrying out the integration over p. The singular properties 
of D^(ki, ei) are similar but related to another pair of spheres: |p + Kki\ — pi = 
and |pi — Kk\ \ —p = 0. The 5-function terms from both of these matrix elements add a 
non-zero contribution to the integral from (fTTj) . As the result, one can express _Rr 12 ] as a 
sum of two terms. The first term, RpL * s only due to the contribution of the 5-singular 
parts of the matrix elements, whereas the second one, -R^] > accoun ts for the integration 
over the whole p-space without the points which belong to the spheres defined above. 
Within the framework of the 'delta'-approximation the term R^ is omitted 



completely, so that instead of i?[i 2 i one uses only R 



[12]" 



Ri 



121 



R 



(S) 
121 



E 



e 2 D 



'l/ 2 U 



k 2 ) +D 



e s + ftwi — E\ 



+ 



+ 



+ 



Ju 2 u 



e 2 



(++), 



e s + hijj x - £] 



s=pi — ?iki 



- S=P2+/UC2 



vpi_ 
2ttH 



2nh 



da 



e 2 D 



U2U 



) (k 2 )/^ ) fn 



p- 



ei) 



e s + fouji - £i 



s=pin p — ftk.i 



41 t (-n p ,e 2 )e 1 D^(k 



C++), 



El 



(12) 



s=p 2 n p +?ik2 

Here the subscript v stands for the set (s, \i) with the momentum s defined as indicated. 
The vector matrix element D^^ a ^(k) (where o"b = '±' and = '±') is given by the 
integral 



Dfc: a) (k) = I dr^t(r) a exp (-ikr) *£">(r) , 



(13) 



and is subject to the conditions |pb + Kk\ 



Kk\ — Ph 7^ 0, which mean 



p a ^ and |p a 

that the matrix elements on the right-hand side of (fT2*j) do not contain the 5-terms. The 
quantities /i^] a (n, e) and b VhV ^(e) are expressed in terms of the unit bispinor u v of a 

[10]: 

t 



plane wave, and the bispinor scattering amplitude Gv ' ( n 



f(±) 



n, e 



(eo)^(n), 



^b^(e) 



ea u L 



(14) 



Accounting for (|T2|) one substitutes the exact amplitude [A^[i2]] , with the quantity 
according to the rule: 



n(J) . • j 
R [l2 ] +1-«[12] , 



(15) 



[M [12] ] + - 

with I [12] defined in (fTUjl. 

Hence, within the framework of the 'delta'-approximation the total 2BrS amplitude 
M. (see ©) acquires the form: 



M 



>(5) 



[21] 



[Ml2]]_+ [M[21)]_ -(16) 
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Here one can account for the approximate formula derived in [5]: 

[M [12] }_+ [M [21] }_ ~ --^| / dr^Ht^e-^+^^+Hr). (17) 

2.3. Application to a point Coulomb field 

In this section we apply the developed approach to construct the approximate amplitude 
of the 2BrS process in the case of a relativistic electron scattering in a point Coulomb 
field, —Ze/r. This case is of interest in connection with the experimental data obtained 
for E\ = 70 keV electrons scattered by various targets [29-32]. It was noted in several 
publications [1,3,4,16,17,41] that for the energies of the electron and the photons within 
the tens of keV range (i.e. as in the experiments) it is the retardation and relativistic 
effects which are mostly important for the formation of the 2BrS spectrum rather than 
the effect of screening due to atomic electrons. 

The additional difficulty (as compared to the non-relativistic case) in applying the 
formalism of the relativistic 'delta'-approximation to the scattering in a point Coulomb 
field appears due to the absence of the closed analytic expression for the scattering 
states wavef unctions. To overcome this difficulty we make another approximation and 
use the Furry-Sommerfeld-Maue wavefunctions, which are accurate up to the order Za 
and for which the closed analytical representation is known (see, e.g. [9,10,42,60]). 
The use of the FSM wavefunctions does not lead to essential difference from the results 
obtained within the framework of the relativistic distorted partial wave approach if the 
main contribution to the cross section comes from the partial waves of high orbital 
momentum I ^> Za [9,42]. 

So far there have been no numerical investigations of the domain of applicability 
of the FSM approximation to the 2BrS process. However, such an analysis was carried 
out for the single-photon BrS process [61] and for the Compton scattering of photons 
from bound electrons [49]. It was established that for high Z elements (Z > 60) the 
calculations based on the FSM approximation underestimate the exact partial wave 
results for double differential cross section by 10-50% while for smaller Z the difference 
is rarely greater than 20% and typically is essentially less. Theoretical estimations of 
the range of validity of the FSM approximation with respect to the atomic number Z in 
connection with construction of the relativistic Green's function can be found in [42,62]. 

The FSM wavefunction is given by (e.g., [9]): 

¥j\r) = e 7 ^ 2 T(l - iaO e ipr / h M - ^ «vj f(kt£, 1, -i(pr - apr)/h} u v . (18) 

Here £ = Zae/pc, the notation F(a, b, z) stands for the confluent hypergeo metric 
function, T(z) is the Gamma-function, u u is the unit bispinor amplitude of a plane wave. 
In ()18|) and in what follows the notations a is used for (a) '±' if a is a superscript, and 
(b) '±1' if it is a factor. 

Using ()18)) in ()13j) one derives the following expression for the one-photon matrix 
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element: 



[eD[X a) (k 



FSM 



The factor Cp^a an d the integrals /, 



PbP. 



PhPa V / 



H ("^(k)) - («[I ( p CT - b) (-k)]*)(ea) 



PbPa 



k), I 



> x PbPa 



(k) are defined as follows 



CjX a) = e^+W/a T(l + ia b £ b ) T(l - ia a £ a ) , 

/KffaJ^j _ y ^ r g-i(Pb-Pa+/ik)r/?i ^p(cr b )* ^(cr a ) ^ 

lK^)(k) = — / dre -i(Pb-Pa+^k)rM F K)* r VF K) 

PbPa V / 2^ / P b \ Pa 



(19) 

(20) 
(21) 

(22) 



Here the short-hand notation Fp° h ^ = F fia^a 6b,a, 1, — i(Pb,a r — 0"b,aPb,a?") /^J is used 



for the hypergeometric function. The three integrals from (|19|) are related through [9]: 

(23) 



where q = pb — p a + Kk. Following the procedure described in [9] , one transforms the 
term l£^^(k) to the Nordsieck-type integral [63] and then derives: 

iC^ftt = _ (g 2 - 2 Pb q + ia b 0)^^ a^-^) 



g 2 (q 2 + 2p a q - ia a 0) i(T ^ 
x jqF(-i(x b £ b , ia a £ a , 1, + ia b £ b F^l - ia b £ b , 1 + ia a £ a , 2, z 



r g 2 (Pb ~ ^bPb n Pa ) - 2q(p b q) 
q 2 - 2p b q 



x 



q2 



where 



? 2 (PbPa - 0"a(TbPaPb) ~ 2(p b q)(p a q) 



(24) 



(25) 



(g 2 -2p b q)(g 2 + 2p a q) 
By analyzing the asymptotic form of the FSM wavefunction (JTSj) one derives the 
following expression for the function /^ a (n, e) (see lfHJ) ): 



(<t) fn el 



/W a (n, e) - a ^ 5(n Pa - an) < (ea) , (26) 

iPa 



where (n, e) is given by: 



i"l (n, e) 



3 i(7| a ln((l-CTCOS a )/2) 



£ a l'( 1 • - irrC) e , , , 

n ; r u' ea 

p a r(l + ia£ a ) l-acos# a Vbl ; 



1 an 

2e a 



Pa 



an 



^a,(27) 



with cos^ a = n Pa n. 

Now we are ready to write the amplitude Ai^ (see 116)1 ) within the framework of 
the FSM approximation. For the negative-energy part one obtains: 

[Mi2]]_+[AV]L eiG2 

where k = ki + k 2 . 



< [^(k) - alfeS(k) - alK^-k)] u vl (28) 
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^ < { w / dfip ^ (k2) - a * (k2) - * ( " k2) ^ 



is aj6j = ab. 

Finally, using lfI5)l and <|2S )l in ([T2 ft . and 



«^(29) 
is adopted 



cap + e p + 7 °m p c 2 ) [a^Cki) - a.al^ki) - alg* (-k x . 
where the subscripts z, j denote the Cartesian coordinates, and the rule which 

carrying out the summation over // one 

+ r 



writes: 

^[2] = 4 e ii e 2i wis, 



^[121 (Pi -^l) 



ft[gj(p 2 + /ik 2 ) 



+ 



ia r(i - ieo r 



+ 









/ sin 8 P 


d9 p / d 




Jo 


J 




r~ 


p2tt 




/ sin 8 P 


d9 p / d 




Jo 


./Q 





7j[jy(p 2 n p + ftk 2 



Here 6* p , P are the poh 

/C [12] 



= 4 a, 



are defined as follows 

^ +) i P ; s +) (k 2 ) + sfc- ) (6)ifcr ) ( k : 
^ +) iS; ) *(-k 2 ) + ^ ) (6)i: 



(30) 

state. The tensors 



ajoc 



i; } *(-k 2 ) 



a«j >A(s, , (31) 



>.i 



= ajA(s, u>i) 

— (XiOL 



OLi 



^ +) 4 ( P1 +) (ki) + < 1 +) (6)4 ( P " 1 +) (k 1 ; 



/7(-+),. 

^spi A spi 



cfc +,i i.t." ) *(- k i) + s«+ +, (&)ii;, +, *(-ki 

ln[(l-cos0i)/2] 



aa 



i c , 



(32) 



= rt~ ) 

P2S 1- COS 01 



a: 



;i / p -- ) (k 2 ) - %aI p ; s -Hk 2 ) - crf^-M^ 



x A(s, ^i)ai 

,i& In[(l-cos6> 2 )/2] 



— S — ^ki' 



= rf(++)_ 

j« spi l-cos# 



1 - — a(p 2 - s + hk 2 

ZE2 



a 



,-A(s,u;i) 



^(kO-al^l-kO^j 

cosf^i = pi(s + Akij/pj and cosy 2 = p 2 (s 
quantities Bplp^\l;) an d A(p,c<j) are given by 

' ' ' 1 ; ' 1 cap + e p + 7°m p c i 



x M+ + Hki) - a.alipT 1 , Fi= , 

Here cos^i = pi(s + fiki)/^ and cos 6*2 = p 2 (s — fik 2 )/p 2 is used in (J3*3~j) and 



(33) 



(34) 
The 



PbPa lU r(l+iC) PbPa ' lP ' ] e p {e p + hj-e l ) 
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Figure 2. Dependence of Z~ 2 A A a ullU2 , f^) on the nucleus charge The incoming 
electron energy, the energies of the photons and the geometry of the emission are 
as indicated. Filled circles with the error bars are the experimental data from [29]. 
Solid lines correspond to the relativistic (thick line) and non-relativistic (thin line) 
'delta'-approximations. Dashed lines are used for the dependences obtained within 
the framework of relativistic (thick line) and non-relativistic (thin line) Born-Elwert 
approximations. Triangles present the results of the exact non-relativistic non-dipole 
approximation [4]. 



with e p = ^Jp 2 c 2 + rripC 4 . 
3. Numerical results 

The approach described in sections 12.21 and 12.31 was applied to calculate the spectral- 
angular distribution of the 2BrS formed in an electron scattering in Coulomb fields of 
a variety of charges. The figures correspond to the incoming electron kinetic energy 
E\ = 70 keV, except for figure El where E\ = 10 keV. In figures El lU and El the energy 
of the first photon is fixed at huj\ = 20 keV, while for the second photon it varies 
within the range hu2 = 2.5. . .47.5 keV. In figures H3 and HI the curves are plotted for 
fojj\ = hjj2 = 25 keV and ftwi = 1 keV, frwi = 1 . . . 8.5 keV, respectively. 
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Nucleus charge Z 



Figure 3. Calculated cross section d 4 <r WlW2 (f2i, O2) as a function of Z versus the 
experimental data from [30] . Solid lines correspond to the relativistic (thick lines) and 
non-relativistic (thin lines) 'delta'-approximations. The results of relativistic plane- 
wave Born and relativistic Born-Elwert approximations are drawn with the short- 
dashed and the long-dashed lines, respectively. The energies £1, Wi ; 2 and the emission 
geometry are as indicated. 

The results in all graphs refer to the com planar geometry when the vectors pi, ki 
and k 2 lay in the same plane. The emission angles (#1,2, ^1,2) , are measured with respect 
Pi. Because of the axial symmetry the differential cross section dV/^doJido^df^df^) = 
dV^o^fii, f2 2 ) depends on the difference (ipi — y^) but not on the angles ipi and y?2 
separately. 

The kinetic energy e± = 70 keV, the photon energies and the geometry of the 
radiation correspond to the experimental conditions [29-32]. 

Numerical calculation of the relativistic amplitudes was carried out within the 
framework of the 'delta'-approximation and using the FSM wavefunctions. The 
analogous non-relativistic calculations were performed with the exact Coulomb 
wavefunctions. The results of relativistic Born calculations presented in the figures were 
obtained by programming the formulae presented in [5] . The curves corresponding to the 
relativistic Born-Elwert approximation represent the relativistic Born curves corrected 
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Figure 4. The calculated cross section d 4 a- WlW2 (rii,r2 2 ) as a function of lo 2 versus 
the experimental data [30] for a 70 keV electron collision with an Au atom (Z = 79). 
Solid lines correspond to the relativistic (thick line) and non-relativistic (thin line) 
'delta'-approximations. Short-dashed line represents the relativistic plane-wave Born 
approximation, the long-dashed one stands for the correction due to the Elwert factor 



by the Elwert factor [60] (we used its relativistic analogue proposed in [64]): 

P! 1 - exp (-27T&) 
p 2 1 - exp (— 27T4 2 ) 

In figures I2HH we compare the results of calculation of 2BrS cross section 
obtained within the framework of relativistic 'delta'-approximation with the available 
experimental data [29-32] and the the results from other theories. The latter 
include the exact non-relativistic non-dipole treatment of the 2BrS process in a point 
Coulomb field [4,17], the non-relativistic non-dipole 'delta'-approximation [3], and the 
relativistic plane- wave Born approximation with and without the correction factor 1)35)1 . 
The calculations within the framework of the non-relativistic dipole approximations, 
although having been carried out, are not presented here. This is because, as it is 
known [3,41], the dipole-photon scheme, which neglects the correction terms of the 
leading order v/c, strongly underestimates the magnitude of the 2BrS cross section. 
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Figure 5. Results from the relativistic 'delta'-approximation theory versus the exact 
non-relativistic calculations [4] for Z~ 2 d 4 cr (i)lW2 (rii, ^2) in a 70 keV electron collision 
with an Ag atom (Z = 47). The sets of curves (the 'delta'-approximation) and symbols 
(the exact non-relativistic results) correspond to different emission geometries: (a) 
solid line and triangles are used for 9\ = 6% = 30°, (b) long-dashed line and squares 
for 9\ = 8 2 = 45°, and (c) short-dashed line and circles for 0x = #2 = 90°. 

One of the conclusions which can be drawn on the basis of the data presented in 
figures EHU is that the relativistic effects result in a decrease of the magnitude of the 
cross section. This is clearly seen if one compares the dependences obtained by using 
the relativistic (thick line) and non-relativistic (thin line) 'delta '-approximations, which 
are presented in figure El Qualitatively, the influence of the relativistic corrections could 
be illustrated in terms of classical electrodynamics as follows. Apart from the spin- 
related effects, the general consequence of a movement with relativistic velocity is the 
dependence of mass of a projectile on its velocity, which leads to the increase of the mass 
asm = m (1 — v 2 /c 2 )~ 1 / 2 (m is the rest mass). This effect reduces the magnitude of 
the projectile's acceleration, which defines the intensity of radiation. 

Presented results show that the relativistic 'delta'-approximation, on the whole, 
reproduces quite well the available experimental data. However, there are two 
exceptions. In figure El in the range Z ~ 100 the calculated values underestimate the 
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Figure 6. Comparison of the results of relativistic 'delta'-approximation with 
the exact non-relativistic calculations with accounting for radiation retardation [17] 
of 2BrS cross section Z~ 2 d 4 <7 Wl(i)2 (Oi, O2 ) as function of energy of second emitted 
photon U2- Presented four graphs correspond to collision of 10 keV electron with four 
different atomic targets with Z = 1, 13, 47, 79. The energy of the first emitted photon 
uj\ = 1 keV and emission geometry 6\ — 62 — 45°, <p\ — f2 = 180° are the same for 
all graphs. Results of relativistic 'delta'-approximation are given by thick solid line. 
For results of non-relativistic exact calculations the thin solid line is used. Relativistic 
plane-wave Born curves are plotted by dashed line. 

experimental data. The possible explanation of this discrepancy is in the inadequacy 
of the FSM approach in the limit aZ ~ 1, The second example of a strong deviation 
between the theory and the experiment is seen in figure H] at the tip-end of the spectrum. 
However, in this case, the discrepancy may be due to the experimental error, since the 
experimental data for U2 = 50 keV lies far above all theoretical predictions. 

In figure we compare the predictions obtained within relativistic 'delta'- 
approximation scheme with the results of the exact non-relativistic theory (the 
retardation included) [4] . It was demonstrated in the cited paper in the non-relativistic 
domain the 'delta'-approximation is in a good agreement with the exact calculations. 
The account for the relativistic effects leads to the decrease in the cross section 
magnitude. This influence has been already mentioned above. We note also that this 
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decrease is dependent on the geometry of the emission. The nature of this feature is 
not fully clear at present, and we believe to get the explanation as soon as the exact 
relativistic results will become available. 

For the lower values of the projectile energy the influence of the relativistic 
corrections is less pronounced. This is demonstrated by figure |B1 where the results 
of the relativistic 'delta'-approximation (thick solid lines) for a 10 keV electron are 
compared with the exact non-relativistic non-dipolar calculations (thin solid lines) 
presented in [17]. Let us point out that the non-relativistic theory based on the 'delta'- 
approximation produces the results (not plotted in the figure) which practically coincide 
with the thin lines for all elements presented. The difference between the thick and thin 
solid curves increases with the atomic number Z. Although the exact nature of the 
difference can be established when the exact relativistic theory of the effect becomes 
available, some part of it can be attributed to the use of the FSM wavefunctions which 
accuracy also decreases with Z. 

For the sake of comparison in figure El we present the curves corresponding to the 
relativistic plane- wave Born approximation (dotted lines). It is clearly seen that the 
simpler theory becomes absolutely inapplicable in the range of medium and large Z 
values. 

On the whole, we may state that the developed approach is adequate for the 
description of 2BrS process and is an efficient tool for numerical analysis of the 
characteristics of the process (the spectral and spectral-angular intensities of the 
radiation). The results obtained within the 'delta'-approximation (both relativistic and 
non-relativistic) describe quite well the behaviour of 2BrS cross sections over wide ranges 
of energies of the projectile electron and the photons and of atomic numbers Z. The 
influence of the relativistic effects increases with e\ and Z, and they must be accounted 
for in addition to the radiation retardation effect. 

4. Conclusions 

Basing on the numerical results presented above, we regard the developed approach as 
more reliable than the simpler relativistic theories. Within the framework of the 'delta'- 
approximation the analytical structure of the 2BrS amplitude is simplified considerably. 
Additionally, this method allowed, for the first time, to carry out numerical analysis of 
the two-photon free-free transitions in relativistic domain beyond the plane wave Born 
approximation. The method is easily generalized to the case of a n-photon free-free 
transition. 

Nevertheless, despite the (relative) simplicity and the efficiency of the outlined 
formalism, by no means does it solve in full the problem of the free-free relativistic 
transitions. The exact relativistic treatment of the process is still unavailable. This 
task is difficult from the analytical and the computational viewpoints. The direct 
evaluation of the amplitude by means of the partial wave expansion of the initial/final 
states wavefunction, the Green's function and the multipole expansions for the emitted 
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photons seems hardly to be implemented in the nearest future. 

However, we are more optimistic on the prospects of two other approaches. The first 
one is based on the exact relativistic treatment of the 2BrS process within the framework 
of the FSM approximation, where one can construct, in addition to the wavef unctions, a 
practically useful form for the Green's function (see [42]). The second approach is based 
on the Sternheimer-Dalgarno-Lewis method [58, 59] for the calculation of two-photon 
transition amplitude between two states of a continuous spectrum. The work in these 
directions is being carried out. 
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